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ABSTRACT 

The emergence of the magnetic field through the photosphere has multiple manifestations 
and sunspot s are the most prominent examples of this. One of the most relevant sunspot 
properties, to study both its structure and evolution, is the sunspot area: either total, umbra or 



penumbra area. Recently Schlichenmaier et al. (2010) studied the evolution of the active region 



(AR) NOAA 11024 concluding that during the penumbra formation the umbra area remains 
constant and that the increase of the total sunspot area is caused exclusively by the penumbra 
growth. In this presentation the Schlichenmaier's conclusion is firstly tested, investigating the 
evolution of four different ARs. Hundreds of Intensitygram images from the Helioseismic and 
Magnetic Imager (HMI) images are used, obtained by the Solar Dynamics Observatory, in order 
to describe the area evolution of the above ARs and estimate the increase and decrease rates 
for umbra and penumbra areas, separately. A simple magnetohydrodynamic model is then 
tentatively used in a first approximation to explain the observed results. 



1. Introduction 

Sunspots are among the most prominent man- 
ifestations of interaction between the solar mag- 
netic field and solar plasma. One of the most rel- 
evant sunspot properties to study is the structure 
and evolution of their area, either umbra, penum- 
bra or total area. 



Hathaway & Choudhary (2008) computed the 



daily values of the total area of one active region 
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during 14 days. Recently, [Schlichenmaier et al. 
f2010) studied the evolution of one active region, 
NOAA 11024. They computed the areas of um- 
bra and penumbra during 6 hours and concluded 
that during the penumbra formation the umbra 
area remains constant and the increase of the total 
sunspot area is caused exclusively by the penum- 
bra growth. 

One of the goals of this work is to compute 
the area evolution of four active regions, during 
their transit through the visible hemisphere, us- 
ing images taken at intervals of approximately 10 
^ 20 min. We then concentrate on the study of 
the area evolution after the active regions have 
reached their maxima, a period that we call de- 
cay phase. In order to explain the observations, 
we propose a simple model whereby a radial fiow 
has a dominant role in the umbra/penumbra evo- 
lution, both through advection and stretching of 
the magnetic field. 



2. Data series and fuzzy area 

In this study we use the images from Solar 
Dynamics Observatory (SDO) with 4096 x 4096 
pixels and time steps of approximately 10 to 20 
minutes. The observed active regions are NOAA 
11117, with observations starting in 26 October 
2011, and NOAA 11428, NOAA 11429, and NOAA 
11430, with observations starting in 8 March 2012. 
Figure [T] shows examples of the kind of images used 
for the study of the considered active regions. 




Fig. 1. — Snapshots of the active regions considered. 
NOAA 11117 in the left panel, NOAA 11428,11429,11430 
in the right panel. "Courtesy of NASA/SDO and the AIA, 
EVE, and HMI science teams." 

Due to the difficulty to define the exact limit 
of the sunspots and in particular to define the 
limit between umbra and penumbra, we use fuzzy 
sets to compute the area of the sunspots. Fig- 
ure |2] shows the intensity histogram obtained from 
a smoothed image, from which the membership 
functions to the fuzzy umbra and penumbra may 
be computed. The degrees of membership of each 
intensity value to the fuzzy umbra and penumbra 
are obtained using values t/i, ^i and P2, which 
delimit respectively the regions of uncertainty cor- 
responding to the separation between the umbra. 



penumbra and photosphere (Fonte & Fernandes 



[2009 J . 
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Fig. 2. — Example of the histogram obtained from one 
smoothed image. Values Ui and U2 represent the limits of 
the region of transition between the umbra and penumbra 
and values Pi and P2 the limits of the transition between 
the p enumbra and the photosphere. ( [Fonte Fernandes | 
[2009t . 



3. Area calculation and simulation 

The areas of the regions corresponding to the 
fuzzy umbra, penumbra and the total sunspot of 
the four active regions indicated below are com- 
puted using a fuzzy approach, which provides in- 
formation on the uncertainty of the obtained val- 



ues (Fonte & Fernandes 112009 



Fonte & Lodwick 



20041). The obtained results are shown in figures 
3) |4j |5] and [6) where the variation of the mini- 
mum and maximum possible values for the area of 
the total sunspot (respectively the red and green 
lines), the umbra (respectively the dark blue and 
purple lines) and the penumbra (respectively the 
light blue and brown lines) are represented. 

Using the information extracted from the evo- 
lution of the active region area, we compute the 
growth rate for the maximum umbra and the max- 
imum penumbra areas and the decay rate for the 
maximum area of the total active region using a 
linear fit. The results are described in tabled) 

In a first approach to interpret the obtained 
results, a simple diffusion- ad vect ion numerical 
model using the MHD formalism was used. In 
this model we simulate a circular formed sunspot 
and extract the area evolution for the decay phase. 

We assume a circular sunspot with only radial 
velocity of the form Ur = u^r/rQ^ where r is the 
radial coordinate and ro is the initial radius of 
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Fig. 3. — Area evolution of the active region NOAA 
11117. 



Fig. 5. — Area evolution of the active region NOAA 
11429. 
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Fig. 4. — Area evolution of the active region NOAA Fig. 6. 

11428. 11430. 



Area evolution of the active region NOAA 



the sunspot. The system of equations to solve in 
cyhndrical coordinates (r^O^z) assuming d/dz = 
d/dO = is represented by equations 1 and 2, 



combination of Bessel functions of the first kind. 
The expressions that we obtain are represented in 
equations |3] and [4j 



dt 



d^B, 1 dB, 



dBr 



^ dr ^ V dr'^ r dr 



dur dBr 
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dr 
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Br 



r dr 



(2) 



where r] is the mag netic diffusivity an d has the 
value of 0.2 km^s-^ ( |Chae et al. 12008 ). 

The initial field considered is constructed based 
on the observation of NOAA 10933 from IBorrerol 



et al. (2011 ), where we adjust the data to a linear 
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where Jo and Ji are Bessel functions of the first 
kind of order zero and one respectively. 

The results obtained in the simulation using uq 
equal to -0.035 Km are represented in table 
[2] and in figure [?) The value for uq is adjusted in 
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Table 1: Results obtained for the growth rate, de- 
cay rate and maximum area for each active region. 



Table 2: Results obtained using the model com- 
pared to the results obtained by fuzzy method us- 
ing the observations. 



NOAA Max area 



{fiHem) 



11117 
11428 
11429 
11430 



Umbra 
growth rate 
{liHem/day) 



Total area 
decay rate 
{fiHem/day) 



800 
480 
1600 
350 



21 
12 

40 
6 



-127 
-85 
-142 
-84 



NOAA 



11117 
11428 
11429 
11430 



Simulated 
decay rate 
{jiHem/day) 



Total area 
decay rate 
{jiHem/day) 



-94 
-61 
-143 
-51 



-127 
-85 
-142 
-84 



order to produce the best agreement between the 
observation and the simulation result. 




Fig. 7. — Example of the results obtained using the 
model. 



4. Conclusions 

Our results seem to be consistent with lSchlichenl - 



maier et al. (2010) results, though not conclusive 



yet. Although the umbra growth rate is much 
weaker than the penumbra's, there is some evolu- 
tion that can be noticed. 

The numerical model tested here indicates that 
the decay rate of the sunspots could be explained 
by advection and stretching of the magnetic field 
due to a negative radial velocity field. 
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